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The U = oo Hubbard model on a hypercubic lattice with exactly one hole is studied. It is seen 
that a degeneracy can occur in the maximal-spin sector. It is proved variationally that the Nagaoka 
(maximal spin) state is not a ground state, if degeneracy exists. We give an explicit construction of 
a variational state, an incommensurate spin- wave state, with energy lower than that of the Nagaoka 
state. For the case of more than one hole, a few examples are discussed where a degeneracy exists. 
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The Hubbard model plays a key role in the study of magnetism, metal-insulator transitions, and high-temperature 
superconductivity. It is a simple model that captures the physics of phenomena of correlated systems. In view of 
the difficulty in establishing rigorous results of the model, the Nagaoka's theorem [jy is very significant, providing 
a good example, though of a restricted use, of metallic ferromagnetism in itinerant electron systems. The result of 
Nagaoka that the ground state of the Hubbard model with an infinite on-site repulsion with exactly one hole is a 
i , maximal-spin state is remarkable, given that for large on-site repulsion with no holes a maximal-spin state cannot be 
a ground state, due to the antiferromagnetic correlations that arise from kinetic exchange phenomenon. There have 
been several contributions extending or simplifying the original result of Nagaoka and investigations on the 

t^- ■ instability of the Nagaoka state, when the number of holes is more than one (5|-[7]] . In the following we give a simple 
proof that the Nagaoka state is not a ground state if there is a degeneracy in the maximal-spin sector ]?J . We follow a 
variational construction of a spin-wave state, and in conjunction with the fact that a degeneracy exists, and show that 
by lowering the spin one can lower the energy below that of the Nagaoka state. The simplest case where a degeneracy 
occurs is the one-hole case for a closed chain with an odd number of sites, and similarly the higher-dimensional lattices 
' with an odd number of sites in at least one direction. 

Let us consider the Hubbard model on a hypercubic lattice with N sites (iVi x N2 x A3.., with A, number of sites in 
■ the i'th direction), and z coordination number, with periodic boundary conditions. The Hamiltonian can be written 
O 1 as 

\o ' 

H=-tr, ]T clc ja +UJ2 mirm (1) 

<ij><7 



where ^(cj^) annihilates (creates) an electron of spin a at i'th lattice site, and = c[ a Ci a is the number operator. 
The first sum in the above equation is for nearest neighbours only, with the hopping strength t. The parameter 77 = 1 
for the case of fermions (that we will discuss through out), and -q = — 1 for hard-core bosons. All the conclusions that 
we draw for the fermion problem carry through to hard-core bosons, by changing the sign of rj. We consider the case 
of the on-site repulsion U to be infinitely large. This forces all the physical states with any number of electrons to 
consist of no doubly-occupied site orbitals, tantamounting to a super-Pauli exclusion. We choose a basis set of states, 
a S^-dimensional Hilbert space, with this constraint built in and consequently the Hamiltonian consists of only the 
first term in the above equation. The total spin S, and the z-component of the total spin S z (at i'th site the spin 
operators are given as S z — (n^ — rijjJ/2, = cj^c^) are good quantum numbers of the model, and we use them 
to label the sectors of the Hilbert space. Let the number of electrons N e = N — 1. The basis states consist of direct 
products of N — 1 site orbitals occupied each with either a spin-up or down electron. Also we note that all the matrix 
elements of the Hamiltonian have a definite sign, by choosing the set of basis states carefully. In the sector with 
S z — Smax — N e /2, the basis states can be labeled by the location of hole, and are obtained by operating with an 
annihilation operator on a reference ferromagnetic state (no holes and no down spins). The basis state i is given by 

\i>=c n \im ~>=c n \REF> (2) 

where the hole is at the i'th site. The reference state \REF >— Yit=i n c l-\\® > ^ s created from the vacuum state 
by using the creation operators. With this choice of the basis set, the matrix elements have all one same sign, 
< i|-ff|j >= tQj where = 1 if the site i and j are nearest neighbours on the lattice, and zero otherwise. For 
the sectors with S z ^ S ma x the construction is similar. For instance, for S z = S max ~ 1 sector, the basis states are 
{\ik >}. A basis state is labeled by the locations of the hole (at i) and the down spin (at k), and is constructed from 
the reference state as \ik >= c^c\^c^\REF >. The overlaps between the basis states are all of the same sign (using 
the fermion anticommutation relations), 

< i'k'\H\ik >=*Cm' +hv$i',k) (3) 
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where the first term in the parenthesis comes from a movement of the up-spin electron, and the second term from 
that of the down-spin electron. It should be noted that the above equation does not change for hard-core bosons (by 
using the boson commutation relations), if we use 77 = — 1 in Eq.(l). 

As can be seen from the above equation, the Hamiltonian in a matrix form looks like H — tA, in any S z sector, 
where A is matrix with positive-semidefinite elements, j4y > 0, and the size of the matrix depends on S z . The 
diagonal elements are zero, and there are exactly z nonzero elements on any row (of magnitude unity). The largest 
eigenvalue is X m ax < z independent of size, from Gershgorin theorem. And its easy to check that a uniform vector 
(which is a maximal-spin state) is eigenstate of the matrix, with eigenvalue z, implying X m ax — z. If the matrix is 
'connected' ^,||]in the configuration space (i.e. the elements of B — A M , for a finite M, are positive definite, Bij > 0), 
Xmax is a non-degenerate eigenvalue for any finite size from Perron-Frobenius theorem. For a bipartite lattice, where 
Ni, N 2 .. are all even, the sign of t is irrelevant. This can be seen by a canonical transform, changing the sign of the 
electron creation and annihilation operators on one of the two sublattices. It suffices to study one sign of t. So for 
t < 0, the lowest eigenvalue of the Hamiltonian corresponds to the highest eigenvalue of the matrix A. In this case, 
the Nagaoka state (S = S max ) is a ground state, with an energy — z\t\, and also it is nondegenerate (connectivity 
of A is provable in more than one dimension) ||). For nonbipartitc lattices, where at least one of Ni,N2.- is an 
odd number, we have to study both signs of t, and the nature of the ground state is quite different, as we shall see 
below. We proceed for t < similarly as above, as the problem boils down to the investigation of X m ax, and the 
Nagaoka state is a ground state. The more interesting case, however, is for t > 0, where a degeneracy exists in the 
maximal-spin sector. Here we are confronted with the lowest eigenvalue X min of A and it very well could depend on 
the size of A; all the sectors of S z have to be examined. Table 1 shows a few examples of periodic two-dimensional 
lattices where a degeneracy occurs, along with the lowest eigenvalues from numerical diagonalizations for two spin 
sectors S = S max and S max — 1, and the Nagoaka state fails to be a ground state. It can be shown from a variational 
construction that degeneracy disfavors high-spin states j?]] . We give a simple proof based on a spin- wave state that 
the Nagaoka state is not a ground state, as a direct consequence of a degeneracy. Our variational state has simple 
structure, and is similar to Roth-type wavefunctions investigated earlier H0. The strategy is to look at two sectors 
with S z = S m ax and S z = S max — 1, and show that X m in{Smax — 1) < X m in(S max ), which implies that a state with 
total spin S = S m ax — 1 has lower energy than the energy of Nagaoka state. 

We focus on a one-dimensional closed chain with an odd-number of sites, and all the arguments used below can be 
trivially extended to higher-dimensional systems. In the sector with S z — S max we have only one species of electrons, 
which makes the system a spinless fermion problem. The term with the coulomb interaction U = 00 in Eq. (1) does 
not come into play at all, and all the many-particle eigenfunctions are constructed as direct products of the single- 
fermion states. The translational invariance of the lattice helps us in diagonalizing the Hamiltonian in momentum k 
space. The energy spectrum (the many-electron spectrum or the one-hole spectrum) is given by Sk — 2t cos k, where 
k = 2im/N and n an integer. For t > (and setting t = 1 from now onwards), the lowest energy of the spectrum is 
for k = —ttzLtt/N with a two-fold degeneracy (and the two degenerate states have a net momentum), and the highest 
energy is for k = which is non-degenerate. Similarly for a square lattice with Ni x N 2 sites and N\ odd and N 2 
even, the lowest-energy state is two-fold degenerate. In this sector all the states carry a total spin of S = S m ax, and 
hence this is an orbital degeneracy. Let us denote the two degenerate lowest-energy states in this sector as 

<pi=^2a t \i>- ) (t) 2 = ^2b l \i> . (4) 

We have H\<j>\^ >= Em\4>i,2 >> where the energy is Em = — 2cos7r/iV, and the following conditions on the eigen- 
functions from normalization and orthogonality 

5^0?* = = l;J2 a *ibi = 0; (5) 

and further since the two states are eigenstates of the Hamiltonian, we have 

<ij> <ij> <ij> 

More importantly the two degenerate eigenfunctions are related by a phase 

aj = bje^ (7) 

where q Q = 2tt/N. As we shall see later this fact will be very crucial and useful. 

We now turn to the sector with S z = S m ax — 1, where an eigenstate of the Hamiltonian will either carry a total 
spin S = S max or S — S max — 1, as no other spin is possible. To diagonalize the Hamiltonian in this sector is a 
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nontrivial problem as we have one down spin this time. However, we may not be impeded by that. The normalized 
eigenstates with S = S max are easy to construct as the Hamiltonian commutes with the total spin, so a lowering of 
the z-component of the states we constructed above gives us the eigenstates we desire here. The two maximal-spin 
lowest-energy degenerate states are S~(f>i^, with an energy Em- The eigenstates with S = S max — 1 are not easy to 
construct. We resort to a variational method. Let us construct a general variational spin-wave state 

Ml) = ^=J2J2a i e i ^\ik> (8) 

where the basis state \ik > has a hole at the site i and a down spin at fc, as discussed earlier. It is clear, from 
Eq.(5), that the above state is normalized. For q — the state is an eigenstate of the Hamiltonian with maximal 
spin. For q 0, the state does not have a definite spin, the expectation value of the spin is less than S max , as it 
has a component from both S max and S max — 1 spin sectors. If q ^ 2irn/N , with an integer n, the above is an 
incommensurate spin-wave state. It is straightforward to evaluate the variational energy of this state using Eqs. (3) 
and (6), and we get 

E a (q) =< 4> a (q)\H\(f> a (q) >= E M + ^ ~ ^ M (9) 
where the fourier-transform function 

A(q) = ]T ataje^*-^. (10) 

<ij> 

It is easy to show that A(q) > Em, by using a variational principle on a state <j) — a i exp(igri)|i > in the S z — S max 
sector, where <f>i (see Eq.(4)) is a ground state. This implies Em < E a {q). The equality holds at two points q = and 
q = q a , as we will see below. From a similar construction of a spin- wave state 4>b(q), using &, instead of a, in Eq.(8), 
we get the variational energy £*(g) = Em + [B(q) — Em]/N — 1, where 

B(q) = £ b*b je ^-^. (11) 

<ij> 

The two fourier-transform functions are related as the two degenerate states are related by a phase (see Eq.(7)); we 
have 



A(q)=B{q-q ). (12) 

The above equation, along with Eqs. (6), (10) and (11), implies A(q ) = B(0) = Em = A(0) = B(—q ). We have 
now two different spin- wave states constructed from the two degenerate lowest-energy states with S = S max . Though 
both these states have energies higher than Em, we can now superpose these states to hope for a lowering in energy, 
as their overlap could be nonzero and useful. A precaution, that we have to take, is that the variational state should 
not become a spin eigenstate with S = S max . Let us try a variational state 

4> = 1(Mq)-MJ)) (is) 

with a choice, from hindsight, q' = q — q Q . Though we are restricting the variational parameter space, we shall see 
below that we still have enough freedom to minimize the variational energy below Em- 

The above variational state has a norm < ip\ip >= N/N — 1, and the overlap between the two components of the 
state for the choice q' = q — q a is < <f)i,(q')\H\(f) a (q) >= B{q) — B(q a ) — Em/N — 1, where we have used Eqs. (3), (7) 
and (11). The variational energy E^ —< il>\H\ip > / < ip\ip > is, using the expressions for E a (q), and Eb(q') from Eq. 
(9), 

- E M + 7^(B(q + qo) + B(q - q ) - 2B(q)) (14) 

where we used Eq.(12) to write all the fourier-transform functions in terms of B(q). Now we seek values of q such 
that 7(g) = B(q + q ) + B(q — q ) — 2B(q) is negative, which implies E^ < Em- We note that the function B{q) is 
periodic (with period 2tt), and is non-constant, B(0) = B(~q ) = Em < 0, and B(ir) = —Em- Then it is guaranteed 
that 7(g) < 0, for some value of g, as the function B(q) necessarily has a local maximum. The expression for 7(g) 
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can be substantially simplified using Eq. (11) as 7(g) = —2(1 — cosq )B(q). Now any value of q for B(q) > would 
suffice to give us the desired result 7(g) < 0. For q = ir, we have 7(71") = —4(1 — cosg Q ) cos (ir/N), and 

2 7T 

E$ = E M + —(l-cosq )cos — <E M . (15) 

This completes the proof that the variational state in Eq.(13), with q = ix,q' = n — q , has lower energy than 
that of the Nagaoka (S = S m ax) state. The fact that we have two degenerate lowest-energy states with maximal 
spin is very crucial. The above arguments can be readily extended for more than one dimension. The structure 
of the proof does not change for higher-dimensional systems. For instance consider a two-dimensional lattice with 
periodic boundary conditions, with N=A 1 x N 2 sites. And let N\ be an even number and N 2 odd. The energy 
spectrum for the hole in S z = S max sector is e(ki,k 2 ) = 2(cosfci + cosk 2 ), with fc, = 27m,- /iVf. The lowest energy 
is for fci = 7r, k 2 — it ± ir/N 2 , with Em = —2(1 + cos(n/N 2 )), with a two-fold degenerate. The various quantities 
q, q ,{ri} now become two-dimensional vectors. The phase between the two degenerate eigenstates is, in Eq.(7), is 
go = (0,2tt/N 2 ) and we choose q = (0,7r) in Eq. (14), to get the desired result < Em- In a general case of a 
hypercubic lattice with N = N± x N 2 x N3... sites, and with at least one of Ni is an odd number, we get a two- fold 
degenerate maximal-spin state for t > 0, where the Nagaoka state is not a ground state. We expect that a similar 
proof can be constructed for any lattice with exactly one hole, when a degeneracy exists. 

We would like to briefly comment on the nature of the actual ground state, if degeneracy exists. We expect 
similar arguments can be constructed to show that the lower the spin the better is the ground state energy, and most 
probably the ground state has the lowest possible spin. A heuristic argument can be given to understand the trend 
of minimal-spin ground state when degeneracy exists. In the non-degenerate case, the ground state eigenfunction 
is positive definite (i.e. a basis can be chosen such that the probability amplitudes of all the basis states in the 
ground state are positive), which is the key point of the Nagaoka state. If we have a degeneracy, the orthogonal 
states cannot be both positive definite, and at least one of them have negative probability amplitudes. A lowering of 
the spin introduces, by itself, negative signs. For instance, any state with S = S z = Smax — 1 is of the form \S = 
S z = S max — 1 >= J2ik Xik\ik >■ It is necessarily accompanied by the constraints (as S + \S — S z — S max — 1 >= 0) 
Sfe Xik = for a fixed i. This implies that the amplitudes come with both signs to satisfy the constraint equations. 
Now it may be expected that these negative signs, and the negative signs as a consequence of degeneracy may conspire 
to cancel out making the eigenfunction as much 'positive definite' as possible, and translate into a kinetic energy 
lowering. And by lowering the spin further, one lowers the energy further. Hence an expected minimal spin ground 
state. 

We discuss the case when the number of holes Nh is more than one. It is easy to check if the maximal-spin sector has 
a degeneracy in the lowest many-hole energy state, as all we need to solve is a spinless fermion problem. We have the 
following interesting and simple examples where degeneracy occurs, (a) Consider the case of a bipartite lattice where 
for Njj — 1 we have no degeneracy and a Nagaoka ground state. For a closed chain with an even number of sites, the 
one-hole energy spectrum in S z = S max sector is e(k) — —2\t\ cosfc. When we have two holes, the ground state in this 
sector is a direct product of single-hole k— space eigenstates (not true for hard-core bosons) with k\ = 0, k 2 = ±2ir/N, 
which is two-fold degenerate. The same holds for Nh even and less than N — 1. And similarly for an Ni x N 2 lattice, 
for many values of Nh the lowest-energy state in the maximal-spin sector is degenerate. In the particular case of 
quarter filling N e = Nh = N/4, for any N, there is a degeneracy, (b) Consider the case of a non-bipartite lattice, with 
degeneracy for Nh — 1 (which implies Nagaoka is not a ground state from the proof outlined above). For example a 
closed chain with an odd number of sites has a degeneracy in the lowest-energy state of the maximal-spin sector when 
Nh is odd and less than N — 1. Similarly for higher-dimensional systems we encounter degeneracy for many values 
of Nh- (c) Consider the case of a long-range hopping model, where the constraint of nearest-neighbour-only hopping 
that we have in Eq.(l) is relaxed, which exhibits macroscopic degeneracy. It is easy to diagonalize the maximal-spin 
sector, and the one-hole spectrum has only two energy levels, with energies —t(N — 1) (non-degenerate), and —t 
(degeneracy N — 1). For t > 0, the lowest-energy many-hole state for Nh < N — 1 is degenerate. For t < 0, there is 
degeneracy for 1 < Nh < N — 1. The singular case of Nh = N — 1 also has a degenerate lowest-energy state, but 
there is only one spin sector S = 1/2. 

In all the examples we considered above there is a degeneracy in the lowest-energy state of maximal-spin sector. We 
expect that the Nagaoka state is not a ground state in these cases. It would be interesting to see if one can generalize 
our simple proof for the above examples 0. Particularly, it should be investigated whether the ingredients that we 
used for Nh — 1, viz. the translational invariance, which is important to diagonalize the Hamiltonian in S z = S max 
sector and get a handle on the degenerate states (Eq.(7)), and the choice of a basis states such that all the overlaps 
have a single sign (Eq.(3)), are essential. Currently an investigation is under progress to extend our procedure in a 
more general context, viz. for non-cubic lattices with more than one hole. 

In conclusion, we have seen that there is a general trend in the nature of the ground state. For Nh = 1, when 
the maximal-spin sector has a non-degenerate lowest-energy state, the ground state is a maximal-spin Nagaoka state, 
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and a degeneracy disfavors higher-spin states and a minimal-spin ground state is expected. An explicit proof, based 
on a variational construction of a spin-wave state, is given to exhibit that the Nagaoka state is not a ground state if 
degeneracy exists on hybercubic lattices, with an odd number of sites in at least one direction. 

We would like to thank K. Hallberg for many useful discussions, and K. Penc for bringing the work of Suto JtJ to 
our attention. 
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Table 1 The lowest eigenvalue E(S) for two spin sectors S = S max , S max — 1 for various sizes for t = 1. 
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